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ABSTRACT 

We present high spatial resolution, medium spectral resolution near-infrared (NIR) 
H- and if -band long-slit spectroscopy for a sample of 29 nearby (z < 0.01) inactive 
spiral galaxies, to study the composition of their NIR stellar populations. These spectra 
contain a wealth of diagnostic stellar absorption lines, e.g. Mgi 1.575 fim, Sil 1.588 fim, 
CO (6-3) 1.619 /im, Mgi 1.711 ahu, Nai 2.207 psm, Cai 2.263 /xm and the 12 CO and 
13 CO bandhcads longward of 2.29 fim. We use NIR absorption features to study the 
stellar population and star formation properties of the spiral galaxies along the Hubble 
sequence, and we produce the first high spatial resolution NIR HK-b&nd template 
spectra for low redshift spiral galaxies along the Hubble sequence. These templates 
will find applications in a variety of galaxy studies. The strength of the absorption lines 
depends on the luminosity and/or temperature of stars and, therefore, spectral indices 
can be used to trace the stellar population of galaxies. The entire sample testifies that 
the evolved red stars completely dominate the NIR spectra, and that the hot young 
star contribution is virtually nonexistent. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: spiral - galaxies: stellar 
content - infrared: galaxies 



1 INTRODUCTION 

Disentangling the stellar populations of galaxies is impor- 
tant to understanding their evolution and the enhance- 
ment of star formation in the Universe, since the inte- 
grated spectrum of galaxies is sensitive to the mass, age, 
metallicity, dust and star formation history of their dom- 
inant stellar populations. Despite this complexity, the ob- 
served optical spectra of galaxies can be accurately repro- 
duced (e.g., Worthey 1994). High-quality observational data 
is paramount to test and improve spectrophotometric mod- 
els, which are used, in turn, to derive the main properties 
of high redshift galaxies. Optical Lick spectral indices have 
usually been used in stellar population studies of local galax- 
ies (Worthey & Ottaviani 1997), but they have problems 
due to blended features and various population of galaxies. 
Near-infrared (NIR, hereafter) spectroscopy has advantages 
compared to optical, because red giant branch (RGB) stars 



* Based on observations collected at the European Southern Ob- 
servatory (ESO), Chile, programs 068.B-0653 and 076.B-0714 
f E-mail: jarkot@utu.fi 



dominate at 2 fim and the wavelength region contains many 
diagnostic stellar absorption lines (e.g., Mgi 1.575 fim, Si 
1.589/xm, CO(6-3) 1.619 fim, Mgi 1.711 fim, Nai 2.207pm, 
Cai 2.263 /im and CO(2-0) bandhead >2.29 fim) which can 
be used as indicators of stellar population in terms of their 
temperature and/or luminosity. The contributions of dwarf, 
giant and supergiant stars can be characterized by different 
depths of their absorption lines and by different shapes of 
their NIR continuum (e.g. Langon et al. 1999). 

Much of the work at NIR wavelengths has focused on 
unusual objects with either active galactic nuclei (AGN) 
or very strong star formation activity. These surveys in- 
clude those of luminous and ultraluminous infrared galax- 
ies (Goldader et al. 1997; Burston, Ward, & Davies 2001; 
Dannerbauer et al. 2005), various types of starbursts (En- 
gelbracht et al. 1998; Vanzi & Rieke 1997; Coziol, Doyon, & 
Demers 2001; Reunanen, Tacconi-Garman & Ivanov 2007), 
Seyfert galaxies (Ivanov et al. 2000; Reunanen, Kotilainen, 
& Prieto 2002,2003; Riffel et al. 2009) LINERs (Larkin et al. 
1998; Alonso-Herrero et al. 2000), and ellipticals (Silva et al. 
2008; Cesetti et al. 2009). However, relatively little NIR 
spectroscopy exists for inactive nearby spiral galaxies. Such 
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surveys are necessary, though, for understanding their star 
formation histories, and for setting a baseline for the NIR 
emission from inactive spiral galaxies. Without such a base- 
line, the contribution of quiescent stellar populations to NIR 
emission in more exotic systems like ultraluminous infrared 
galaxies and starbursts will remain unknown. 

The two largest previous NIR spectroscopic studies of 
inactive spirals were performed by Mannucci et al. (2001) 
and Bendo & Joseph (2004). Mannucci et al. (2001) stud- 
ied 28 galaxies with a 7 x 53 arcsec aperture and R = 400 
- 500. Generally, they found a high degree of uniformity 
for the spectra of spirals (earlier types being more homo- 
geneous) but their study was hampered by the low spectral 
resolution which prevented them from detecting faint and/or 
blended lines. More recently, Bendo & Joseph (2004) stud- 
ied 41 galaxies in the K-band and 20 galaxies in the H-band, 
with slightly higher spectral resolution R = 800-1200. They 
extracted spectra from 30x15 arcsec apertures, and found 
that in almost all cases evolved red stars completely domi- 
nate the nuclear stellar populations, while young stars were 
virtually non-existent. We stress, however, that both these 
studies suffer from poor spatial resolution, leading to a pos- 
sible mix of bulge and disk stellar populations. 

To improve this situation, we present in this work NIR 
spectroscopy of nearby inactive spiral galaxies along the 
Hubble sequence, to study their stellar population content 
and SF histories. Stellar emission and absorption features 
can probe with great detail the age and SF properties of the 
stellar population in the nucleus and as a function of radius. 
Comparison of these properties along the Hubble sequence of 
spirals will assess any relationship of SF with Hubble stage. 

We present NIR spectroscopy for a sample of 29 nearby 
(z < 0.01) spirals with reasonable spatial resolution (~1 
kpc/arcsec) and medium spectral resolution (R ~600) to 
study the NIR properties and ages of their stellar popula- 
tions. Except for a few individual sources, these are the first 
medium resolution NIR spectra of a sizable sample of low 
redshift spirals. Long-slit spectra of the spirals are extracted 
across the nucleus to detect all the important diagnostic ab- 
sorption lines, and emission from the galaxies can usually 
be studied to large distance from the nucleus. Our aim is 
to study the stellar population of spirals based on the NIR 
stellar absorption indices, and to inter-compare these indices 
within the different morphology types. 

This paper is organized as follows: In Section[2]we de- 
scribe the sample selection, observations, data reduction and 
methods of analysis. In Section^ we present the results and 
discussion concerning the properties of the galaxies, and the 
composite quiescent spectra, and in Sectionf3]we summarize 
our conclusions based on the full sample of 29 spiral galaxies. 
Throughout this paper, we adopt a concordance cosmology 
with H = 70 km s" 1 Mpc" 1 , fi m = 0.3 and fi A = 0.7. 



2 SAMPLE SELECTION, OBSERVATIONS, 
DATA REDUCTION AND ANALYSIS 

This project started as a followup to our NIR spectroscopy 
of Seyfert galaxies (Reunanen et al. 2002, 2003). We created 
a comparison sample of non-active spiral galaxies that were 
matched in Hubble type with our sample of Seyferts. We 
selected in each Hubble type from E to Sc a similar number 



of nearby bright galaxies (z < 0.01; m(K) < 11 to keep the 
exposure time reasonable to reach the required S/N) and 
paid careful attention that the spiral galaxies had no indi- 
cation of nuclear or peculiar activity, based on their multi- 
wavelength properties, especially from optical spectroscopy, 
radio and X-ray data. We planned to observe five galaxies 
for each Hubble type between E and Sc, in order to ob- 
tain a good estimate of the intrinsic spread in the spectra of 
each class, to remove possible individual peculiar character- 
istics. Although our sample is not complete, we consider it 
as representative of the class of nearby spiral galaxies, and 
an improvement over the samples in Mannucci et al. (2001) 
and Bendo & Joseph (2004). 

The H- and K-band spectroscopic observations of the 
spirals were carried out in February 2002 and March 2006 
at the 3.5 m ESO New Technology Telescope (NTT) using 
the 1024x1024 px SOFI camera and spectrograph (Moor- 
wood et al. 1998) with pixel scale 0.288 arcsec pix -1 . The 
average seeing during the observations was ~0.7 " FWHM. 
The medium resolution (R ~600) long-slit HK-bsmd spectra 
were taken along the major axis of the galaxies and across 
the nucleus, with the red grism and slit width 1.0 arcsec. 
The wavelength range of the spectra is from 1.5 to 2.5 pm 
and the useful slit length ~2'. The spectra of each target 
were taken with pairs using ABBA observing cycle of alter- 
nating galaxy and sky exposures nodding ~ 60 arcsec along 
the slit, and slightly moving the target position in the slit 
between exposures, with total integration time of 48 minutes 
(except 32 minutes for NGC 5612), divided into individual 
exposures of 60 s. 

A literature search was carried out to collect previous 
measurements of the galaxies in our sample: intrinsic ve- 
locity dispersion a, and optical spectral indices Mg2 and 
H/3 in the Lick system (e.g., Targer et al. 1998 and ref- 
erences therein). The velocity dispersion and the Mg2 are 
taken from the HyperLEDA compilatiorQ (Paturel et al. 
2003), which has applied unification of the measurements 
from different sources (i.e. Golev & Prugniel 1998). An av- 
erage H/3 was taken if more than one measurement was avail- 
able. The measurements are listed in Table[T] together with 
other properties of the observed sample of galaxies (redshift, 
K-band magnitude, morphological type, Hubble type), and 
the FWHM seeing during the observations. The literature 
references are given in Tabled 

Data reduction was performed using IRAF0 and fol- 
lowed the procedure described in Hyvonen et al. (2009). 
The AB pairs of spectra were subtracted from each other 
to eliminate background sky emission and then divided by 
the flat-field images, constructed from incandescent lamp 
spectra. Bad pixels and cosmic rays were masked out. One- 
dimensional spectra were extracted from each individual im- 
age. Wavelength calibration was done using Xe-arc lamp cal- 
ibration frames using the same continuum tracing as for the 
galaxy, and its zero point was checked against the OH sky 
lines. The galaxy spectra were divided by a G-type spectro- 
scopic standard star and then averaged. The spectroscopic 

1 http : //leda ■ univ-lyonl . frj] 

2 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 1. General properties and journal of observations of the spirals. References for the literature measurements of <r, Mg2, and H/3 
are listed in Table [2] 
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Columns: (1) target name; (2) redshift from NED; (3) the ii"s-band magnitude from 2MASS in a 10 arcsec diameter aperture; (4) and 
(5) morphology and type from NED and RC3, respectively; (6) Hubble type from RC3; (7) velocity dispersion; (8) and (9) Lick Mg2 
and H/3 index, respectively; (10) observing date; and (11) average seeing during the observations. 



standards were selected to be as close to the target galax- 
ies as possible. The nuclear ('nuc', 1 arcsec diameter, corre- 
sponding to ~1 kpc at the average distance of the galaxies) 
and off-nuclear extended emission spectra ('offl' and 'off2', 
from radii 0.5 to 2 arcsec at opposite sides of the nucleus, 
and covering an area of 1.5 arcsec 2 ) were extracted for each 
target. Three or even more bin widths were added together 
to trace the extended emission. The //-band spectra were 
normalized by fitting the continuum level to both sides of 
the absorption features. In the A"-band, the continuum was 
fitted only to the shortward wavelengths from the CO band- 
head. The journal of observations of the sample is presented 
in Tabled] 

The equivalent widths (EW) of the spectral lines were 
measured on the normalized spectra, integrating over the 
same passbands for all targets; A=l. 5870-1. 5910 [im for Sil, 
A=l. 6175-1.6220 ^m for CO(6-3), A=2. 2040-2. 2107 ^m for 
Nal, A=2. 2577-2. 2692 ^m for Cal, and A=2. 2924-2. 2977 )im 
for CO (2-0). The continuum at the line was defined by a 
straight line fitted to the continuum around the spectral 
feature. For purpose of a direct comparison, the same wave- 
length ranges were adopted as in Origlia et al. (1993) and 
Silva et al. (2008). 



The accuracy of the EWs is limited by the uncertain 
continuum level in the vicinity of the particular line rather 
than the S/N of the spectra. This is especially true for the 
//-band which is rich in spectral features, and for the CO(2- 
0) bandhead in the A"-band because the continuum level has 
to be extrapolated longward of 2.295 fim. It is worth to keep 
in mind, that there are many blends and broad molecular 
features (Wallace & Hinkle 1996). Furthermore, the mea- 
sured EWs of the indices are evidently upper limits due to 
spreading caused by velocity dispersion. The possible depen- 
dence of the EWs on velocity dispersion and resolution was 
checked by degrading the medium resolution stellar spec- 
tra from Wallace & Hinkle (1997) to the lower resolution 
applicable for the spiral galaxy sample, and measuring the 
indices of all the five diagnostic features for both resolutions. 
This comparison showed that the indices do not significantly 
depend on the degradation: all measurements agree within 
the errors. Therefore, the bandpasses are wide enough to 
fully encompass the measured features, despite the broad- 
ening due to the intrinsic velocity dispersion. We estimate 
our typical EW errors to be: ~0.5 A for the //-band and 
for the /f-band metallic lines, and ~1.0 1 for the /f-band 
CO (2-0) feature. 
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Table 2. References for the central velocity dispersion, Mg2 and 
H/3 measurements from the literature. Only galaxies from our 
sample with available measurements are listed. For convenience, 
the reference codes are identical with the ones used in Hyper- 
LEDA (Paturel et al. 2003). 
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NGC 5134 


IFC96 


IFC96 
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References: 51057 - Wegner et al. (2003); 7Sam - Faber et al. 
(1989); Be02 - Beuing et al. (2002); BG091 - Bettoni et al. (1991); 
CDB-93 - Carollo et al. (1993); DBD93 - Do Souza et al. (1993); 
DS83 - Dressier & Sandage (1983); ENEARc7 - Bernardi et al. 
(2002); HSMP99 - Heraudeau et al. (1999); IFC96 - Idiart et al. 
(1996); JMA96a - Jablonka et al. (1996); LCO-LO - Davies et al. 
(1987); OFJSB - Dalle Ore et al. (1991); S83 - Schechter (1983); 
Ton84c - Tonry (1984); Tr98a - Trager et al. (1998); VegaOl - 
Vega Beltran et al. (2001); WFG92a - Worthey et al. (1992); WK 
- Whitmore & Kirshner (1981); WM - Whitmore &; Malumuth 
(1984). 



3 RESULTS AND DISCUSSION 

The inactive spiral galaxies have strikingly similar NIR nu- 
clear spectra as a function of their Hubble type, especially 
when compared to more active classes of objects. The main 
H- and K-band spectral absorption features from the stel- 
lar population - Mgl, Sil, CO(6-3), Nal, Cal, and CO(2-0) 
- are clearly visible. Their relative strengths are quite sim- 
ilar for all Hubble types. No H-band [Fe II] or K-band Br7 
or H2 emission lines are evident at a significant level. How- 
ever, weak H2 and/or Br7 emission was detectable in a few 
galaxies after subtracting an elliptical galaxy template rep- 
resenting the old stellar population, (see Section [331> . 

The H K-band nuclear spectra for the individual spiral 
galaxies in the sample are presented in Fig. [1] The EWs of 
all the diagnostic absorption lines were measured for all the 
spiral galaxies, and they are presented in Tables [3] and [4] 
The off-nuclear and nuclear EWs of the spirals are generally 
similar, and they are in broad agreement with the EW values 
from the literature (Mannucci et al. 2001; Bendo & Joseph 



2004; Silva et al. 2008; Cesetti et al. 2009; Hyvonen et al. 
2009), if differences in aperture sizes and index definitions 
are taken into account. 

3.1 NIR line indices as stellar population 
indicators 

The NIR spectroscopy is well-suited to study stellar pop- 
ulations in galaxies and star clusters because it contains 
many strong stellar absorption features tracing the pres- 
ence of red stars, particularly red supergiants and giants. 
The 7f-band has been widely studied (e.g., Merrill & Ridg- 
way 1979; Kleinmann & Hall 1986; Wallace & Hinkle 1997, 
Ivanov et al. 2000), and it is a natural choice for studying 
reddened (dusty) galaxies at low redshift (e.g., Mobasher 
& James 2000; Mannucci et al. 2001). However, the dust 
emission that tends to dilute the stellar absorption has even 
smaller effect in the //-band which has been studied system- 
atically only recently (e.g., Origlia et al. 1993; Oliva et al. 
1995; Ivanov et al. 2004; Fremaux et al. 2007). The J-band 
remains relatively unexplored in galaxies (but see Wallace 
et al. 2000 for stars). 

The integrated NIR light is dominated by different stars 
from those that contribute significantly at visible wave- 
lengths: Whereas the integrated light at visible wavelengths 
comes from a mix of stellar types, the JT-band light in all 
but the youngest stellar systems is dominated by evolved gi- 
ants and asymptotic giant branch (AGB) stars, whose pho- 
tometric and spectroscopic properties are sensitive to age 
and metallicity (Maraston 2005). 

The relative strengths of the NIR features can be used 
to characterize the stellar contents of the central regions 
of galaxies, and to provide insights into their history (e.g., 
Vazquez et al. 2003) because of their sensitivity to the spec- 
tral type (temperature) and/or luminosity (gravity) of the 
star (Kleinmann & Hall 1986; Ali et al., 1995; Ramirez et al. 
1997; Ivanov et al. 2004). Ivanov (2001) discussed the pos- 
sibility to develop NIR analog of the optical H/3 vs. Mg2 
diagnostic diagram used to derive ages and metallicities. To 
give an example for the utility of this diagnostic tool: while 
H/3 vs. Mg2 can be used to constrain the properties of the 
host galaxies of optical supernovae (e.g, Hamuy et al. 2000), 
NIR diagnostic is needed to do the same for the host galaxies 
of obscured supernovae (e.g, Mannucci et al. 2003). Ivanov 
(2001) pointed that the optical H/3 has an IR counterpart - 
Br7 at 2.167 /im, and that there are a number of candidates 
to replace the Mg2, including some NIR Mg features. 

The most evident absorption features in the K-band are 
the molecular CO(2-0) bandhead (>2.295/xm), atomic Br7 
(2.166 /im), Nal (2.207 /im), and Cal (2.263 /tm) lines. Br 7 
is the strongest absorption feature within stars earlier than 
K5, but disappears in stars later than K5 (Kleinmann & 
Hall 1986), while Nal is the strongest atomic feature when 
effective temperature is T cff <3400K (Ali et al. 1995). The 
_ff-band has a very complex absorption line pattern due 
to a number of metallic and molecular lines, especially in 
cool late-type stars. The most studied H-band features are 
atomic Sil (1.589 /im) and molecular CO(6-3) (1.619 /im) 
(e.g., Origlia et al. 1993). The CO is very strong in young 
giant and supergiant stars (10 Myr-100 Myr) and strong in 
cool AGB stars (100 Myr-1 Gyr; Oliva et al. 1995), while 
it is weaker in older population. This makes it suitable to 
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Figure 1. The HK -band nuclear spectra of the individual spiral galaxies. 



trace recent star formation in galaxies (Mayya 1997). There 
are also some relatively weak Mgl features (e.g. 1.504, 1.710 
jj,m) that potentially could provide a direct comparison with 
the optical a element abundances, but they are rather de- 
manding in terms of the quality of the observations (Cesetti 
et al. 2009). The Nal, Cal and CO features in K and M stars 
become stronger with lower temperature, i.e. redder J—K 
colour (Ramirez et al. 1997; Forster Schreiber 2000; Frogel 
et al. 2001). At a given J—K, giants in more metal rich clus- 
ters have stronger Nal and Cal features (Frogel et al. 2001). 
Some line ratios, such as EW(CO 1.62)/EW(Sil 1.59) and 
EW(CO 1.62)/EW(CO 2.29), are even better temperature 
indicators than single lines (Origlia et al. 1993; Ivanov et al. 
2004), especially if the two lines have very similar wave- 
lengths making their ratio insensitive to dilution from dust 
or hot stars. Last but not least, most of these features suf- 



fer some degree of contamination from other elements, much 
like in the optical (e.g., Worthey et al. 1994). 

Currently, there are few self-consistent theoretical spec- 
tral synthesis models for the interpretation of NIR spectra 
of galaxies (although see e.g. Leitherer et al. 1999; Bendo & 
Joseph 2004; Maraston 2005; Riffel et al. 2008). Thus, one 
has to rely on e.g. their comparison with high-resolution NIR 
stellar spectral atlases (e.g., Wallace & Hinkle 1996) or on 
direct comparison with optical indices (e.g. Silva et al. 2008, 
Cesetti et al. 2009). 



3.2 NIR-to-optical and NIR-to-NIR line index 
relations of spiral galaxies 

The optical Mg2, Fe, and H/3 indices show tight correla- 
tions with the central velocity dispersion a (e.g. Bernardi 
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Figure 1. continued. The H if -band nuclear spectra of the individual spiral galaxies. 



et al. 1998; Mehlert et al. 2003), suggesting that the chemi- 
cal and dynamical evolution of galaxies are intertwined. For 
the galaxies in our sample, the relation with a is well estab- 
lished for Mg2 but it suffers from somewhat more scatter for 
the Fe5270. None of the NIR features show as clear a rela- 
tion with a as the Mg2. The best Mg2 analog in the NIR 
appears to be Nal at 2.207 (im which shows approximately 
the same scatter with a as the Fe5270. We also find that 
the Nal at 2.207 fim feature shows a fair direct correlation 
with Mg2 despite the possible combined production of Na in 
both low and high mass stars, which sets it apart from Mg, 
which is produced only in massive stars. 



3.3 Dilution of the NIR line indices of spiral 
galaxies 

Non-stellar thermal sources, such as dust surrounding young 
star forming regions, and/or non-thermal dilution reduce the 
intrinsic values of EWs originating from the stellar popula- 
tion. However, the dust observed in NIR is close to its sub- 
limation temperature, and in such cases HII emission line 
series are expected to be observed. The comparison of ab- 
sorption indices in different filters depends on the level of 
dilution by the non-stellar component that can be differ- 
ent in H- and if-band. The dilution fraction of the contin- 
uum emission at 1.6 /im and 2.3 /xm can be estimated from 
the plots of EW(CO 1.62) vs. EW(CO 1.62)/EW(Si 1.59), 
and EW(CO 1.62) vs. EW(CO 1.62)/EW(CO 2.29), shown 
in Fig. [2l top panels. Objects with diluted stellar features 
lie away from the loci occupied by stars, along the dilution 
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Table 3. EWs of the H-band absorption lines for the nuclear and off-nuclear spectra of the spiral galaxies in units of A. 
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Column (1) gives the name of the spiral galaxy; (2) the luminosity class T; and (3) to (14) the EWs of the main diagnostic absorption 
lines in units of A, in the nuclear 1.0 arcsec diameter aperture (nuc), and in the two off-nuclear apertures (offl and off2). 



vectors, and the fraction of non-stellar continuum is simply 
given by the displacement of the point in the diagram (e.g., 
for dilution from non-thermal emission in AGN see Fig. 3 
in Ivanov et al. 2000). Fig. [2] shows that in the H- and K- 
band, all spirals have either no or only slight dilution. This 
indicates that no significant amount of dust is present in the 
nuclear region of spirals. 

Fig-El bottom left panel, shows the EW(Nal+Cal) as 
a function of EW(CO 2.29) for the nuclear spirals, and for 
elliptical galaxies (from Hyvonen et al. 2009), together with 
stars of different spectral type. Caution should be exercised 
when comparing stellar and galactic data, as galaxies are 
composite stellar systems and their spectra show the inte- 



grated contributions from stars spanning a range of prop- 
erties. Still, given that the NIR spectral region is domi- 
nated by evolved (RGB) stars, this issue is not as criti- 
cal as at visible wavelengths, where stars contribute over 
a much larger range of evolutionary states (both RGB and 
main sequence stars (e.g., Maraston 2005). All three lumi- 
nosity classes of stars have a relatively tight linear correla- 
tion between EW(Nal+Cal) and EW(CO 2.29). Most of the 
off-nuclear spirals are in good agreement with the relation 
of giant (III) stars, while ellipticals deviate more strongly 
from that relation. As expected, for spirals in most cases 
the main role is played by class III giant stars, possibly with 
some contribution from class V dwarf stars. 
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Table 4. EWs of the K-band absorption and emission lines for the nuclear and off-nuclear spectra of the spiral galaxies in units of A. 
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Column (1) gives the name of the spiral galaxy; (2) the luminosity class T; and (3) to (23) the EWs of the main diagnostic emission 
and absorption lines in units of A, in the nuclear 1.0 arcsec diameter aperture (nuc), and in the two off-nuclear apertures (offl and off2). 



On the other hand, the plot of EW(CO 2.29) vs. 
EW(CO 1.62)/EW(Si 1.59) (Fig.d bottom right panel) in- 
dicates that the stellar populations are mostly dominated 
by a mixture of giants (III) and supergiants (I). However, 
it is very likely that the CO (1.62) line is underestimated 
and/or the Sil line overestimated decreasing the line ratio 
of EW(Sil)/EW(CO 1.62) and thus, in both cases, shifting 
the real data points to the left in Fig. [2l bottom right panel. 
The Sil line is slightly overestimated, because it is contami- 
nated, mainly by CO, due to the intrinsic velocity dispersion 
of the galaxies, together with the medium-resolution of our 
spectra. In conclusion, the majority of spirals are closer to 



giants and supergiants than the ellipticals, demonstrating 
the difference in their stellar populations. 

Almost all spirals and ellipticals have EW(CO 2.29) > 
15 A which is larger than the EW of K giants (from 10 
A to 15 A) or main sequence stars. In fact, the only stars 
with EW of CO(2.29) higher than 15 A are M giants and K 
and M supergiants. The stellar population of these galaxies 
must therefore be dominated by very cool stars, probably 
M giants, to match the measured EWs. Note, however, that 
since M stars are much brighter in the K-band than K stars, 
this does not imply that M giants are the only significant 
luminosity-weighted population. 

The EWs of CO 2.29 ^m and Cai of the spiral and 
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Figure 2. For the nuclear spectra of the spiral galaxies (this work) and elliptical galaxies (from Hyvonen et al. 2009): Top left: EW(CO 
1.62) vs. log[EW(CO 1.62)/(EW(Si 1.59); Top right: EW(CO 1.62) vs. log[EW(CO 1.62)/(EW(CO 2.29); Bottom left: EW(Nai 
2.207 + Cal 2.263) vs. EW(CO 2.29); and Bottom right: EW(CO 2.29) vs. log[EW(CO 1.62)/(EW(Si 1.59). The individual galaxies 
are marked according to their Hubble type T (from -5 to 4; E to Sc). Different luminosity types of stars are shown as small symbols: 
supergiants I (open circles), giants III (asterisks) from Ramirez et al. (1997) and dwarfs V (open squares) from Ali et al. (1995). The 
lines in the first two panels enclose the area occupied by stars with no dilution and the arrow gives the direction of the effects of dilution. 
In the last two panels, the lines show the relations for different luminosity types of stars. 



elliptical galaxies are plotted against EW(Nal ) in Fig. [3] 
Majority of the spirals and ellipticals occupy a locus which 
is consistent with the relation of purely old galaxies (Silva 
et al. 2008). They also have larger Nal indices than those 
of solar metallicity cluster stars, consistent with early-type 
galaxies in the Fornax cluster (Silva et al. 2008). Unlike the 
sample of Cesetti et al. (2009), our spirals and ellipticals 
do not show a significant correlation between EW(Cal) and 



EW(Nal), probably due to the smaller abundance range. We 
also note that Ca and Na are produced by both high and 
low mass stars (e.g. Denisenkov & Ivanov 1987; Denissenkov 
2005; Langer et al. 1993; Cavallo et al. 1996), so that yields 
and chemical enrichment history may be relevant. 

Different star formation histories within a galaxy may 
set up a radial stellar population gradient, which can be 
seen as a difference in line strengths and EWs of absorption 
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Figure 3. Top: EW(CO 2.29) vs. EW(Nal). Bottom: EW(Cal) 
vs. EW(Nal). The symbols are the same as in Fig. [2] The dotted 
and long-dashed lines represent the relations of the solar metallic- 
ity cluster star fits and purely old population (Silva et al. 2008), 
respectively. 



lines between nuclear and off-nuclear spectra. In non-active 
galaxies, EWs usually increase toward the bulge where they 
remain roughly constant. In the sample considered here, only 
a few galaxies show any indication of systematic radial vari- 
ation in the strength of the lines between off-nuclear and 
nuclear spectra. While these could be due to a stellar popu- 
lation gradient, deeper and higher spatial resolution spectra 
are needed to confirm the possible gradient. 



3.4 Template spectra 

The spectra of the inactive spiral galaxies indicate that the 
NIR emission from these objects is dominated by old red 
stars that dominate the central bulges of these galaxies. 



By combining each individual redshift-corrected spectra to- 
gether, we can create composite quiescent spectra for each 
Hubble type. These composite spectra can then be treated as 
templates representing an old red stellar component in the 
spectra of galaxies where either continuum emission from an 
AGN, hot dust, emission lines from young stars, or super- 
novae are present. These spectra therefore find applications 
as template spectra for the underlying galaxy which can 
be subtracted from the spectra of composite stellar systems 
(e.g. AGN host galaxies). 

To construct the quiescent template spectra, we divided 
the sample into six groups according to their Hubble type. 
We first redshift-corrected all the spectra to the rest frame. 
Next, we resampled the spectra using splines, normalized the 
spectra and averaged them using the same method used to 
normalize and combine the spectra for individual galaxies. 
The resulting composite quiescent spectra, with the most 
prominent metal absorption lines and the CO band heads 
identified, are presented in Fig. [4] where we also show the 
sky emission and the sky transmission spectra, to give the 
reader an idea of the reliability of the template spectra. 

The differences between the spectra of the galaxies of 
the same class are due both to observational uncertainties 
and to intrinsic differences between the observed galaxies. 
Mannucci et al. (2001) and Bendo & Joseph (2004) found 
inactive galaxies to show a high degree of uniformity, which 
we confirm with our higher spatial resolution data. 

Fig. [5] shows the average EWs of the main absorption 
lines as a function of the Hubble type T. We find a slight 
trend for the EWs of Sii 1.589 pm and CO 1.619 pm to in- 
crease, and the EW of Nal 2.207 /im to decrease, with in- 
creasing T, but the trends are not statistically significant. 



3.5 Spiral galaxies with nuclear and off-nuclear 
emission lines 

The nuclear NIR spectra of three galaxies (all Sb spirals) in 
our sample (NGC 2566, NGC 2935 and NGC 5612) exhibit 
the H2 2.122 /j,m emission line. NGC 2935 is an interesting 
case with red nuclear colour, indicating possible nuclear ac- 
tivity. However, the strength of its absorption lines closely 
matches the quiescent composite spectra derived above. This 
implies that any recent star formation activity (or obscured 
AGN activity) lies deeply buried within an older population 
of evolved stars. 

Weak emission from Br7, the strongest hydrogen recom- 
bination line in the HK-bwads, was detectable in eight sam- 
ple galaxies. Although ionized gas may be present, the neg- 
ligible recombination line flux in the majority of the sample 
implies relatively few young stars, and therefore relatively 
little recent star formation activity. Of the eight galaxies 
with detectable Br7, five are classified in RC3 as Sb spirals 
or later types (T ^ 2), and none are earlier than late SO 
galaxies (T ^ -2). This agrees well with earlier results pre- 
sented in Mannucci et al. (2001) and Bendo & Joseph (2004), 
and supports the paradigm that the relative strength of star 
formation increases along the Hubble sequence, probably re- 
flecting the change along the Hubble sequence in the relative 
dominance of bulge stars to the H- and 7f-band emission. 

Galaxies IC 4214, NGC 2196, NGC 2775, NGC 3715, 
NGC 3892 and NGC 5806 exhibit Br 7 line but nondetectable 
H2 line. The recombination line indicates that young ioniz- 
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Figure 4. The //-band (left-hand panels) and the X-band (right-hand panels) composite nuclear spectra for the spiral galaxies, ordered 
by decreasing Hubble type from top to bottom. The Hubble type is listed in the X-band panels. The locations of the most prominent 
diagnostic spectral absorption features together with the expected locations of main emission lines ([Fell], H2 and Br7) are marked. 
The panels in the third row from bottom show all the composite spectra grouped together, colour coded as in Fig[2] to highlight their 
similarity across the Hubble type. The two bottom panels show the sky emission and the sky transmission spectra, respectively. 



ing stars are present, but the lack of shock excitation lines 
suggests that the stellar population has not evolved to a 
point where supernovae are generating shocks in the inter- 
stellar medium. Therefore, using Starburst99 models (Lei- 
therer et al. 1999), the stellar population must be less than 
3.5 Myr old. On the other hand, NGC 5612 exhibits H 2 
line from shock excitation but very weak or nonexistent Br7 
emission. The absence of strong recombination lines indi- 
cates that they have at least evolved past 8 Myr. If the 



shock excitation lines are interpreted as originating from 
supernovae, then we conclude that the stellar populations 
in these galaxies are in the 8-36 Myr age range. 

4 CONCLUSIONS 

We have presented medium resolution NIR HK-band spec- 
tra of a sizable sample of low redshift inactive spiral galax- 
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Figure 5. The average EWs of the main absorption lines as a 
function of the Hubble type T, for the nuclear spectra (solid lines), 
and for the off-nuclcar spectra (offl and off2, long-dash and dot- 
dash lines, respectively). The average errors in the average EWs 
are ~0.5 A for the metal lines and the CO(1.62), and ~1.0 A for 
the CO(2.29). 



ies, to study their stellar populations along the Hubble se- 
quence, based on the diagnostic stellar absorption lines. Our 
main results are: (i) the EWs of the Sil, CO(1.62) Nal, Cal 
and CO(2.29) features in the spirals are most consistent 
with those of giant stars, whereas ellipticals (from litera- 
ture) show a contribution from main sequence dwarf stars, 
and on average, EW(C02.29) of spirals is somewhat greater 
than that of ellipticals. Most likely, the EWs are not signif- 
icantly affected by non-stellar dilution, (ii) Overall, evolved 
red stars completely dominate the NIR spectra, while con- 
tribution from hot young stars is virtually nonexistent, (iii) 
We produce high spatial resolution NIR HK-band template 
spectra for spirals along the Hubble sequence, that will form 
a baseline for comparison to galaxies with more exotic star 
formation or AON activity. 
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